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We report a numerical study of the tunnel conductance through the Schottky barrier at the contact 
between a semiconducting carbon nanotube and a metal electrode. In a planar gate model the 
asymmetry between the p-doped and the n-doped region is shown to depend mainly on the difference 
between the electrode Fermi level and the band gap of carbon nanotubes. We quantitatively show 
how the gate/nanotube distance is important to get large on-off ratios. We explain the bend of 
the current versus gate voltage as the transition from a thermal-activation region to a tunneling 
region. A good agreement is obtained with experimental results for carbon nanotubes field-effect 
transistors. 



Recently, we have reported logic circuits with car- 
bon nanotube transistors.EJ Single-wall carbon nanotubes 
(CNs) are novel quantum wires consisting of tubular 
graphite sheets p which can be synthesizecLui structures 
~ 1 nm in diameter and microns long.cfu The one- 
dimensional (Id) nature of CNs requires careful analysis 
with respect to the screening and the band bending near 
the contact to the electrode. The purpose of this study 
is to calculate the transmission probability through the 
Schottky barrier at the interface between a semiconduct- 
ing CN and a metal electrode and to compare it with the 
experiments. 

The contact between metal electcad^ a*i^£^|Jjas. 
been studied theoretically before.ffl#»BByfciMN 
The transmission at .-a-xontact has been evaluated by 
ab mifi^jCalculationtlijQ as well as by using simple 
models.Btlu3 Using the latter, it has been shown that the 
charge transfer in CNs is dramatically different than in 
2d or 3d systems and that Ihis-strongly modifies the band 
bending near the contact £30 The band bending varies 
slowly with a logarithmic dependence on the distance in 
a coaxial gate model whepfirthe gate electrode is a cylin- 
der surrounding the CN.I130 The variation in distance 
depends on the gate voltage, but is on the order of the 
radius of the cylinder. For the same reason, p-n junc- 
tions of Id CNs have been shown, to behave differently 
than conventional p-n junctions. E3 

Experiments on samples with semiconducting tubes 
acting as -field-effect transistors have been reported 
previouslyrj'H The Schottky barrier was studied by 
transportE3 and scanning probe experimental^. Typi- 
cally, an oxidized Si wafer was used as the gate. The dis- 
tance between the gate and the CN typically was long, on 
the order of the CN length. Transconductance measure- 
ments showed that the tubes were p-doped, and that the 
capacitive coupling between the CN and the gate was too 
weak to n-dope the CN. Recently, various groups have 
improved the sample device and have reported strong 
doping from_n- to p-type by using a v£rv_close gate of 
oxidized ALEJiarge diameter nanotubes,E3E!j or annealed 
nanotubes Ejt3 

In the present work, we report calculation of the con- 
ductance as a function of the gate voltage in a semicon- 




ducting nanotube. We compare this with experimental 
data. The sample layout is characterized by a planar gate 
that lies very close to the CN. The gate thus strongly 
screens the Coulomb interaction between the electrons 
in the CN. We compare the measurements with a semi- 
classical model based on Poisson's equation. Our method 
of calculation is similar to that in ref.Eil, but our model 
incorporates a planar gate and a CN connecting ste 
to a bulk, electrode, instead of the coaxial gate,l!3! 
electrodeQ, or weakly contacted CN at the interface. 
By using this realistic model, we evaluate the current ver- 
sus gate voltage and compare that quantitatively to the 
real device. A good match is obtained. The numerical 
results of our calculation are shown in detail. Finally, we 
discuss how we can improve the device in order to get 
better transistors characteristics by changing the mate- 
rial and the geometry of the device. 

In our devices, the gate consists of a microfabricated 
Al wire with a well-insulating native AI2O3 layer, which 
lies beneath a semiconducting CN that is electrically con- 
tacted to two Au electrodes on CN. In this configuration, 
the AI2O3 layer of a few nm thickness is much shorter 
than the separation between the contact electrodes (~ 
100 nm). Capacitance measurements on two large Al 
films separated by the same aluminum oxide layer gives 
a thickness of 2 nm, whereas ellipsometry measurements 
give a value of around 5 nm. Figures |l| (a), (b) and (c) 
show measurements of the current as a function of the 
gate voltage for three different samples. The bias voltage 
is V = 5 mV and the measurements are done in vacuum 
and at room temperature. Strong doping from p to n is 
achieved for all samples. Interestingly, an asymmetry is 
observed. There are some sample dependences: We can 
see a smaller on-off ratio of 10 4 in Figs. || (b) and (c), 
compared to the large on-off ratio 10 5 in Fig. fil (a). The 
asymmetry of Fig. [j] (a) is stronger than Fig. jl] (c). The 
current in the valence band is gradually decreasing and 
then rapidly decreasing in Fig. m (a) , but the decrease is 
more slowly in Figs, [j] (b) and (c). We will fit these ex- 
perimental data by numerical calculations based on our 
model and explain how these sample dependence can be 
understood from our model. 

We assume an end-bonded contact model in the cal- 



2 



io- 9 



t 10 10 r 
ZJ 

o 



10 12 



(a) 




-2 2 

Gate Voltage V g [V] 



10-1 2 



f=0.007 
AW/e(0)=0.9 
D =0 
v\ R S =78R 



-4-2 2 

Gate Voltage V„ [V] 




; y 1 

ii 




(d) 








ol.n: 


CN 




■ ^ 










''MB* 



-4-2 2 

Gate Voltage V g [V] 

FIG. 1: (a)-(c) Experimental I(V g ) data for three typical 
nanotube transistors. The experimental results are shown 
by dashed lines. The results from the calculation are shown 
by solid lines, (d) Model of the semi-infinite planar gate at 
z — and y > 0. The source ( or drain ) electrode is an 
infinite plane at y = 0. 



culation, which means that the electrons tunnel in the 
end of the tube from the metal contact. Indeed, Bock- 
rath et al. have shown that the CN is cut into segments 
when the electrode is patterned on top of the CN and 
that transport involves tunneling into the ends of the 
nanotubeo We used a planar-gate model as shown in 
Fig. [I] (d) . A CN is surrounded by dielectric material of 
dielectric constant k = 5, with a distance R s from the 
gate. The Poisson equation relates the potential 4>(y) at 
the surface of CN to the Id charge density p{y) in the 
CN and the potential ^>(y) of the gate, 
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with the Fourier transformation (j) q and of the poten- 
tials 4>(y) and ^(y), respectively, and, 

U g = - {I (qR)K (qR) - K (2qR s )} , (4) 

K 

M q = cxp[-\q\R s ], (5) 

where Iq and Kq are the modified Bessel functions. The 
first term of the kernel Eq. (|J) describes the self- 
capacitance within the tube. The second term is the 
term of mutual capacitance, i.e., the coulomb potential 
due to the charge on the gate, which is induced by the 
charge on the CN. The analytic expression is a good ap- 
proximation for the case R s 3> R with a radius R. In 
calculation of the mutual term, we approximately pre- 
sume the CN to be a wire, instead of taking into account 



the dependence on z of the charge density and the poten- 
tial 4> on the surface of CN. Furthermore, the interface 
between the insulator and vacuum in the experiment is 
ignored in the model. We can neglect the correction of 
an image charge in the insulator at the opposite side of 
CN in the second term of Eq. (||), because it is shown 
in the calculation that the second term is much smaller 
than the first term. 

In equilibrium, the charge density is related to the en- 
ergy E (y) — E n (y) — Ep of the charge neutrality level 
Eq of the CN measured from the Fermi level Ep. We 
obtain 



p(y) = po(y) 



8irR 
Via 2 



D a E () ( y )e-y/ d + ^\ sign(y), 



(6) 

with the lattice constant a, the pinning strength Dq, an 
effective decay constant d of the surface states that equals 
about 2 nmjlj the doping fraction / (the number of doped 
charge carriers per atom of the CN)tL3, and 

Po(y) = e J dev(e)sign(e)F({e - £ (y)}sign(e)) , (7) 

where F(e) — 1 / (exp(e / fceT) + 1) is the Fermi distribu- 
tion function. Equation (j^) is valid if Eq (y) varies slowly 
on the scale of the Fermi wavelength. 
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Tie density of states v in semiconducting CNs is given 
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with the step function 9, the discretized wave number 



Kn =R-( n -l 



(9) 



the band, parameter 7 = -\/3(i7o/2, and the transfer inte- 
gral 70 H3 

Due to the conservation of the total electron energy 
the charge neutrality level is related to the electrostatic 
potential Eq. (|l]) , 



E (y) + e4>(y) = AW, 



(10) 



where AW is a constant and one of the fitting param- 
eters. Note that, although AW relates to the work 
function difference between bias electrodes and CNs, it 
strongly depends on the surface and randomness near the 
contact. 

We solve Eqs. (Q), (0), and (^) self-consistently to 
obtain the potential (fyy) with a boundary condition 
4>(0) = 0. In order to ensure the boundary condition, 
we use antisymmetric sources 



= V g siga(y), 



(11) 



and A VFsign(y) instead of the right hand side of Eq. (|T(]) . 
For this purpose, image charges at the opposite side of 



3 



the interface between CN and the metallic electrode have 
been added in Eq. (||). 

Transmission coefficients for a single barrier are cal- 
culated in the WKB approximation. The transmission 
coefficient T n to the states of the n-th band in the CN is 
given by 



T n (E) = cxp 




E-Epjy) 

6(0) 



dy 



(12) 

where the half of the band gap e(0) = 7/3i?, the integral 
is taken in the barrier, and the energy E is measured from 
the Fermi energy. The calculations show that the con- 
tribution of the higher sub-bands is negligible compared 
with that of the lowest band, because of the effectively 
huge barriers for these sub-bands. For small bias volt- 
age V, the current is given by the Landauer formula at 
finite temperature by the use of the calculated potential. 
A correction for impurity scattering is ignored, since the 
transmission probability in the CN is much larger than 
the tunneling probability at the Schottky barrier. We 
can ignore interference within CNs at room temperature, 
and the use Ohm's law for the series resistance of the two 
barriers near source and drain electrodes. 

We now calculate the transistor characteristics for the 
CN devices. In the following numerical calculations, a 
number of parameters are fixed at the experimental val- 
ues. The energy gap is fixed at 2e(0) = 0.7 eV, whichia 
the average value measured on nanotubes with STM.EZl 
By Joking the experimentally deduced value 70 = 2.6 
eV,E£l this results in a radius R = 0.53 nm. The 
calculations are done for a nanotube with a length of 
A = 380i? = 200 nm and surrounded in a uniform di- 
electric with k = 5. The bias voltage is set at V = 5 
mV and the temperature at T = 300 K. Note that the 
horizontal position of the I(V g ) curve is a free parameter 
in our fit. We make this choice because experimentally 
the whole I(V g ) curve can be shifted horizontally in a 
hysteresis way when voltages as large as 4V are applied 
on the gate. This hysteresis is also observed in samples 
with a Si gate and possibly originates from trapped oxide 
charges within the insulators. 

Figures |l] (a), (b) and (c) show that the calculated 
I{V g ) are in good agreement with the measurements. The 
model reproduces the gap as well as the asymmetry on 
the p and the n doping regions. In the fit of our model 
we have used four parameters: R s , AW, / and D . R s 
and AW are the most important parameters. The gap 
depends principally on R s and the asymmetry on AW. 
The variation of the parameters / and Dq is used for the 
fine tuning. The parameter / has an influence on the 
overall slope of the curves, and Dq on the slope far from 
the gap. The values of the fitting parameters are given 
in the inset of Fig. [I]. 

Figure || shows the band bending as a function of the 
position from the electrode to the middle of the CN for 
different gate voltages. The bands have been calculated 
with the parameters used for the fit in Fig. [I] (a). The 





Ev : 


1 1 


L (a) 

1 


) 

Position 


100 

(units of R) 




i 

l i 




-(c) 


Ev 


I 


Position 


100 

(units of R) 





1 




e c : 


E F 






e v : 


L (b) 


■ 



100 
Position (units of R) 




100 
Position (units of R) 

FIG. 2: Calculated bottom of the conduction band and top of 
the valence band from the contact to the middle of CN for gate 
voltages V g = -2.3 (a), -0.5 (b), 1.2 (ch and 2.1 (d). The 
corresponding current is shown in Fig. H of AW = 0.65e(0). 
The arrow in (c) illustrate the effect of a higher subband. 



different band diagrams show that V g successively shifts 
the Fermi energy from the valence band to the gap, and 
then to the conduction band. It is thus possible to elec- 
trostatically change the doping of the nanotube over the 
full range from p to n doping. 

The Schottky barrier height is not expected to be 
pinned at the middle of the gap, as it is usually the case 
in 3d, and is that expected for an unpinned junction. 
The reason is that, as has been shown .theoretically, the 
effect of the metal-induced gap statesEa is limited, be- 
cause they are on the ring at the contact, instead of a 
plane for a conventional bulk junction.E^I Here, the bar- 
rier height is equal to e(0) — AW for p-doped CNs [Fig. 
I (a)], and e(0) + AW for n-doped [Fig. | (c) and (d)]. 
The barrier is lower and thinner in Fig. |3 (a) than (c) 
and (d), because the Fermi level of the electrode is close 
to the top of valence band. This gives the asymmetry in 
the I(V g ) curves between the two regions. 

Figure |^ shows also that the band bending varies very 
smoothly in space. The strongest band bending varia- 
tion occurs at a length scale of the order of the separa- 
tion between the CN and the gate, in agreement with 
the ca lcula tions where the gate is modelled by a coaxial 
metabll^O Superposed on this variation, weak kinks can 
be observed. These kinks are situated at the integer mul- 
tiple of e(0); an example is indicated by the arrow in Fig. 
U (c). After the next sub-band is occupied, the screening 
becomes stronger. The kinks are due to the divergence 
of the density of state at the sub-band edges. 

Fig. H shows the variation of I(V g ) for different AW 
values. The parameters R Sl f and Dq are taken equal to 
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FIG. 3: AW dependence of the current for R s /R — 48 and 
Do = 0.1. The experiment shown in |l| (a) can be fitted by 
the calculation for AW = 0.65e(0). 



FIG. 4: Calculated current for several gate distance R 3 , and 
AW = 0.65e(0), Do = 0.1. The dashed line on the left shows 
the thermal-activation slope of — e/fcaT. 



what was obtained in Fig. [I] (a). The arrows (a),(b),(c) 
and (d) indicate the gate voltages for which the band 
diagram was represented in Fig. ^j. The current is sur- 
prisingly larger for situations (b) than for (c) and (d). 
Although (b) corresponds to the case that the Fermi level 
lies in the gap, a non-zero current passes through the tube 
due to the broadening of the Fermi distribution function 
at 300 K. In (c), the Fermi level lies in the conduction 
band, but the large barrier hardly allows for electrons to 
tunnel. The current therefore is lower than 10 -12 A, the 
lower limit in the measurements. 

The calculated current is also shown for other AW. 
The on-off ratio increases with AW, because the Schot- 
tky barrier is thin and low for large AW. This indicates 
that the Fermi energy should be located as close as possi- 
ble to the valence band edge in order to get larger on-off 
ratios. The maximum current and the asymmetry also 
increase with AW. If AW is negative, the calculated cur- 
rent is the same after the mirror inversion of gate voltage 
around V„ = -0.68 V. 

Figure H shows the dependence of the I(V g ) curve on 
the gate distance R s . These calculations show the im- 
portance of the gate-CN separation for the on-off ratio. 
For thinner insulator thickness, a larger on-off ratio is 
obtained. Indeed, a thin insulator thickness gives a thin 
barrier width and thus a larger current. For the same 
reason the gap increases with R s . The current corre- 
sponding to the conduction band is out of range in Fig. 
I for R s = 200R. A very large gate voltage V g = 13.4 V 
is needed to obtain a current of 10 -12 A. In previous ex- 
periments, the separation between the CN and the gate 
was large. It indeed did not allow to obs erve a current 
corresponding in the n-doping region.l 16 !! 17 ! 

Further calculations (not presented here) show that 



the slopes around V g = —0.5 V depend on temperature, 
but poorly on parameters like R s , AW or Do (for small 
enough R s , say R s < 100i?). The temperature depen- 
dence in this V g region originates from the broadening 
of the Fermi distribution function at room temperature. 
The current can be approximated in a thermally acti- 
vated form / oc exp[— eVg/ksT] with Boltzmann constant 
kg, as shown in the figure by a dashed line. For small 
R s values, the band bending is mainly located near the 
interface. This changes for large R s . Then, the band 
bending variation becomes comparable to the nanotube 
length and the exponential dependence changes. 

The separation R s between the CN and the gate as 
found in the fit appears to be larger than 48R. This 
value is in relatively poor agreement with the experimen- 
tal value that is of order 10i?. This discrepancy probably 
originates in the description of the bias electrodes. In the 
experiment, the electrodes may better screen the electric 
fields from the gate to the CN than in the model. The 
assumption that the electrons tunnel in the end of the 
CN may also be too strong. The electrons are in fact 
likely to tunnel into the CN at a point that is not right 
at the electrode edge, but close to it. 

In conclusion, we have studied the tunnel conductance 
through the Schottky barrier at the metal-CNs interface 
with a planar gate model. Good agreement with the ex- 
periment is obtained for several samples. The asymmetry 
between the p-doped and the n-doped region has been 
shown to depend mainly on the difference between the 
electrode Fermi level and the band gap of carbon nan- 
otubes. We have shown that the band bending variation 
in the nanotube length is long. We have also shown that 
the choice of source and drain electrodes is important in 
order to obtain a large on-off ratio. We have explained 
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the bend of the current versus gate voltage as the tran- 
sition from a thermal-activation region to a tunneling 
region. 
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